UNCLASSIFIED 

ao  431844 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawingB,  speci¬ 
fications  or  other  data  are  used  for  smy  purpose 
other  than  in  connection  vlth  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


CATALOGED  BY  DDC 


00 

rH 

00 


General  Research 

Confined  Vortex  Flows  with  Boundary-Lajcr  Interaction 


20  FEBRUARY  1964 


Prepared  by  M.  V.  ROSENZREIG.  W.  S.  LEAVELLEN,  and  0.  H.  ROSS 
Aerodynamics  and  Propulsion  Research  Laboratory 


Prepared  for  VICE  PRESIDEN  T  and  GENERAL  MAN  AGER 

1  AHORATORA  OPERATIONS 

00 

196 * 

T*  -1 

M  ROSRU’I  C 

CO 

C-  a(  -  ro  . 


REPORT  NO. 

ATN-64  (9227)  2 


)R  I’OR  \  I  I  'IN 


Report  No. 
ATN-64(9227)-2 


GENERAL  RESEARCH 


Confined  Vortex  Flows  With  Boundary- Layer  Interaction 


Prepared  by 

M.  L.  Rosenzweig,  W.  S.  Lewellen,  and  D.  H.  Ross 
Aerodynamics  and  Propulsion  Research  Laboratory 


AEROSPACE  CORPORATION 
El  Segundo,  California 


20  February  1964 


Prepared  for 

VICE  PRESIDENT  AND  GENERAL  MANAGER 
LABORATORY  OPERATIONS 


( 


ABSTRACT 

Axisymmetric  flow  of  an  incompressible  fluid  in  a 
right-cylindrical  vortex  tube  bounded  by  planar  end 
walls  is  considered.  The  mutual  interaction  of  the 
primary  vortex  flow  field  and  the  end-wall  boundary 
layers  is  studied  by  relating  the  stream  function  and 
circulation  in  these  regions.  An  interaction  parameter 
is  defined  which  determines  the  magnitude  of  the  in¬ 
fluence  of  the  boundary  layers  on  the  circulation  and 
mass-flow  distributions  in  the  primary  flow.  These 
distributions  are  obtained  numerically,  primarily  as 
functions  of  the  interaction  parameter  and  the  Reynolds 
number  based  onthe  radial  flow.  Experimental  results 
are  interpreted  in  the  light  of  this  theory,  and  it  is 
concluded  that  substantial  turbulence  levels  (though 
somewhat  less  than  estimates  made  without  consider¬ 
ation  of  boundary  layer  effects)  must  be  assumed  to 
explain  the  experimental  results. 
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NOMENCLATURE 
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A 

c 

F 


P 

Q 

°BL 
Re 


Re 


Re. 


r 


r 

6 


Boundary-layer  interaction  parameter  defined  in  Eq.  (17) 
Coefficient  of  friction 

Dimensionless  part  of  the  stream  function  which  is 
independent  of  z 

Dimensionless  part  of  the  stream  function  which  is  linear 
with  z 

Length  of  the  vortex  chamber 
Pressure 

Volume  flow  per  unit  length 

Volume  flow  in  the  end-wall  boundary  layer 

Radial  Reynolds  number  =  pur/p 

"Effective"  turbulent  radial  Reynolds  number 

Tangential  Reynolds  number  =  pvr/p. 

Radial  coordinate 

Radial  velocity,  positive  outward 

Tangential  velocity 

Axial  velocity 

Axial  coordinate 

Velocity-profile  shape  parameter  in  the  boundary  layer 
Circulation  -  vr 

End-wall  boundary-layer  thickness 
Dimensionless  parameter  defined  in  Eq.  (5) 
Dimensionless  radial  coordinate  squared 
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NOMENCLATURE 


(Continued) 


P 

v  = 

£ 

p  = 

4, 

Subsc  ripts 
0 

e  = 

m  = 


Viscosity 

Kinematic  Viscosity 
Dimensionless  axial  coordinate 
Density 

Stream  function  defined  in  Eq.  (7) 

Denotes  value  at  outer  edge  of  the  vortex  tube 
Denotes  value  at  the  edge  of  the  exhaust  hole 
Denotes  minimum  value 
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I.  INTRODUCTION 


Vortex  tubes  of  one  form  or  another  have  been  of  interest  to  fluid 

dynamicists  for  sonic  time.  Applications  such  as  the  Ranque-Hilsch  tube, 

cyclone  separator,  and,  more  recently,  the  MHD  vortex  power  generator* 

and  gaseous-core  nuclear  rocket,^  have  provided  the  stimulus  for  research 

in  this  field.  Analytical  and  experimental  investigations  of  confined  vortex 

flows  have  been  numerous.  The  theoretical  treatments  have  been  concerned 

with  the  main  body  of  rotating  fluid  and  have,  for  the  most  part,  been 

restricted  to  one -dimensional  flow  models  (i.e.,  tangential  and  radial 

velocities  being  functions  of  radius  only).  Donaldson^  examined  the  full 

Navier-Stokes  equations  under  these  assumptions  and  found  a  family  of  exact 

solutions.  These  exact  solutions  are  not  readily  applicable  to  the  flow  in  an 

actual  vortex  tube,  however,  because  of  their  inability  to  satisfy  boundary 
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conditions  imposed  by  the  end  walls  and  exhaust-hole  geometry.  Others  ’ 
have  developed  approximate  theories  in  which  the  radial  distributions  of  axial 
velocity  are  specified  (rather  arbitrarily)  in  order  to  satisfy  conditions 
imposed  by  the  exhaust  ports.  A  notable  exception  is  a  study  by  Anderson^ 
in  which  the  laminar  end  wall  boundary-layer  interaction  problem  was  first 
treated. 

Lewellen  has  reviewed  these  studies  and  has  clarified  the  mathematical 
foundation  of  the  approximate  solutions.  In  particular  he  has  shown  that  under 
certain  circumstances  (which  are  usually  met  experimentally),  the  one- 
dimensional  hypothesis  can  yield  a  valid  first  approximation  to  the  solution. 

In  regions  where  this  is  the  case,  the  axial  velocity  distribution  is  determined 
to  the  same  degree  of  exactness  by  its  boundary  values.  In  an  actual  vortex 
tube,  these  are  determined  by  the  boundary  layers  on  the  end  walls  and  by  the 
location  and  geometry  of  the  exit  port  (or  ports),  as  mentioned  above.  Solu¬ 
tions  based  on  arbitrarily  specified  axial  velocity  distributions  have  been 
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reasonably  successful  in  predicting  the  qualitative  behavior  of  the  measured 

tangential  velocity  distributions ,  numerical  discrepancies  being  attributed 
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usually  to  turbulence  in  the  flow  field.  ’  ’ 

The  boundary  layers  which  are  formed  on  the  end  walls  of  the  vortex 
chamber  have  been  studied  previously  without  regard  to  their  effect  on  the 
primary  vortex  flow  which  drives  them.  Boundary-layer  growth  on  a  finite 
disc  under  the  influence  of  a  vortex-like  outer  flow  has  been  studied,  origi¬ 
nally  by  Taylor  ,^and  more  recently  by  Mack^  ,  King,  ^  Weber,  ^  Rott,  ^ 
and  others.  These  studies,  in  which  the  external  velocity  distributions  were 
specified  functions  of  the  radial  coordinate,  have  demonstrated  that  such 
boundary  layers  are  capable  of  transporting  large  quantities  of  fluid  radially 
inward.  This  radial  flow  is  coupled  through  the  continuity  equation  to  an 
axial  flow,  so  that  in  fact  it  is  the  boundary-layer  flow  which  determines, 
for  the  most  part,  the  axial  velocity  distribution  in  a  vortex  tube. 

For  a  fixed  total  flow  through  the  tube,  diversion  of  fluid  into  the 
boundary  layers  reduces  the  radial  flow  through  the  main  body  of  the  vortex 
and  tends  to  weaken  the  vortex  strength,  or  circulation.  In  an  actual  vortex 
tube  (such  as  pictured  in  Fig.  1),  however,  the  boundary  layers  must  at  some 
point  discharge  their  contents.  If  this  process  is  accomplished  in  such  a  way 
as  to  return  the  boundary-layer  fluid  to  the  main  flow  before  it  has  lost  much 
of  its  angular  momentum  and  before  it  is  passed  axially  out  of  the  tube,  then 
(as  is  described  in  this  paper)  it  is  possible  to  maintain  a  near-constant 
circulation  distribution  in  the  vortex.  This  point  is  illustrated  in  Fig.  2: 
fluid  following  paths  1  and  2  will  contribute  to  the  support  of  the  circulation 
distribution,  while  that  following  paths  3  and  4  will  not. 

Mass  flow  may  be  caused  to  leave  the  boundary  layer  for  two  reasons. 
First,  the  bounda ry- layer  theory  itself  indicates  mass  ejection  when  the 
external  circulation  distribution  decreases  with  decreasing  radius.  ^ 

Second,  certain  discontinuities  in  the  wall  geometry,  such  as  a  circular  step 
or  the  sharp  edge  of  an  exhaust  orifice  (i.e.  ,  regions  where  boundary-layer 
theory  is  not  locally  applicable),  can  also  induce  mass  ejection.  Experimental 
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evidence  of  such  boundary-layer  mass  ejection  is  plentiful;^’  ^  several 
examples  are  given  in  Figs.  3  and  4.  Figure  3  shows  boundary-layer  mass 
ejection  in  the  vicinity  of  the  exhaust-hole  radius,  resulting,  most  likely, 
from  both  mechanisms.  Figpre  4  illustrates  the  effect  of  a  circular  step  in 
the  end  wall  at  approximately  mid-radius.  Mass  ejection  is  again  observed. 

The  full  three-dimensional  flow  pattern  in  a  vortex  tube  is  thus  the 
result  of  a  complicated  interaction  between  the  primary  vortex  flow  and  the 
end-wall  boundary  layers.  The  purpose  of  the  present  paper  is  to  study  this 
interaction  theoretically  in  order  to  demonstrate  the  unusual  mutual  influences 
which  are  exerted  by  the  primary  and  secondary  flows  in  right-cylindrical 
vortex  tubes  of  the  type  illustrated  in  Fig.  1.  The  results  are  compared  with 
experiments  and,  as  a  corollary,  ce rtain  observations  are  made  concerning 
turbulence  levels  in  confined  vortex  flows. 
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II.  DESCRIPTION  OF  FLOW  MODEL 


An  incompressible,  constant  viscosity  fluid  is  assumed.  For  the 
purposes  of  the  analysis  the  flow  is  considered  divisible  into  three  regions, 
as  indicated  in  Fig.  1.  It  will  be  assumed  that  the  ratio  of  volume  flow  to 
circulation  is  small  so  that  in  Region  I,  the  primary  flow  region,  the  tangen¬ 
tial  velocity  is  a  function  of  radius  only.  Region  II  consists  of  the  boundary 
layers  on  both  top  and  bottom  end  walls,  in  which  axial  velocity  gradients 
and  therefore  axial  viscous  stresses  become  significant.  Region  III,  bounded 
by  an  imaginary  cylinder  of  radius  equal  to  that  of  the  exhaust  hole,  includes 
the  flow  in  the  vicinity  of  the  axis  of  the  tube,  where  both  tangential  and 
radial  velocities  vanish. 

Appropriate  solutions  are  employed  in  each  of  these  regions  and  a 
consistent  flow  picture  is  obtained  by  application  of  suitable  matching  condi- 
tions.  The  solution  of  Lewellen,  '  which  determines  the  circulation  distribu¬ 
tion  in  terms  of  the  stream  function,  is  used  in  Regions  I  and  III.  This 
solution  becomes  exact  in  the  limit  of  large  circulation  and  requires  that  the 

stream  function  be  specified  at  the  periphery  of  the  applicable  regions.  The 

14 

approximate  boundary-layer  solution  of  Rott  is  used  along  the  end  walls, 
up  to  the  edge  of  Region  III,  to  ,.  'vide  a  second  relationship  between  the 
stream  function  and  circulation  in  the  primary  flow,  thereby  establishing  the 
interaction.  The  mathematical  framework  is  completed  with  the  specifica¬ 
tion  of  the  stream  function  in  Region  III.  Provision  is  included  to  allow  the 
consideration  of  boundary-layer  mass  ejection  due  to  the  exhaust-hole 
discontinuity.  Since  little  of  a  quantitative  nature  is  known  about  this  process 
at  the  present  time,  this  phenomenon  is  treated  parametrically.  The  fraction 
£  of  boundary-layer  mass  flow  at  the  radius  of  the  exhaust  which  is  returned 
to  the  main  flow  before  leaving  the  tube  is  introduced  into  the  computational 
program  and  allowed  to  vary  from  zero  to  one.  The  complete  set  of  equations 
is  solved  iteratively  using  a  digital  computer  for  a  range  of  the  governing 
parameters  . 
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Before  proceeding  with  the  analysis,  several  comments  are  in  order. 
First,  the  division  of  the  flow  into  regions  as  indicated  is  not  equivalent  to  the 
usual  division  of  flows  into  a  viscous  boundary  layer  and  an  inviscid  outer 
flow.  In  the  present  problem,  the  entire  flow  field  is  considered  viscous; 
however,  in  the  boundary  layers  the  axial  velocity  gradients  predominate, 
while  in  Regions  I  and  III  only  the  radial  gradients  are  considered. 

Second,  the  details  of  the  flow  in  Region  III  are  oversimplified  as  a 
result  of  insufficient  knowledge.  Recirculation  of  fluid  originating  downstream 
of  the  exhaust  hole,  which  frequently  occurs  along  the  axis  of  the  vortex  tube, 
has  been  neglected.  This  can  only  be  justified  by  the  belief  that  such  a 
phenomenon  does  not  strongly  influence  the  velocity  distributions  in  Regions 
I  and  II. 

Finally,  a  real,  and  possibly  important,  phenomenon  is  not  considered 
at  all  in  this  paper,  i.e.  ,  the  boundary  layer  on  the  cylindrical  containing 
walls  and  the  details  of  the  jet-mixing  region  if  the  vortex  is  driven  by  tangen¬ 
tial  fluid  injection,  as  is  commonly  the  case.  These  effects  may  influence  the 
end-wall  boundary-layer  growth  and  will  certainly  make  the  determination  of 
the  peripheral  circulation  difficult  in  an  actual  experiment.  For  the  present 
purposes,  however,  it  is  assumed  that  the  outer  cylinder  is  rotating  and 
porous  so  that  such  questions  may  be  avoided. 
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III.  DERIVATION  OF  THE  INTERACTION  EQUATIONS 


Consider  a  cylindrical  coordinate  system,  r,  $,  z,  with  associated 
velocity  components  u,  v,  w.  The  continuity  and  Navier-Stokes  equations 
in  these  coordinates  for  an  incompressible  fluid  with  constant  viscosity  are: 

9(  ru)  .  9(rw)  =  0  (1) 

9  r  <Tz 


9u 

u— +  w 


du 

0z 


2 

v 

r 


_I®E  + 

p  9r 


_1  fhj 
r  9r 


(2) 


u 


9  v 

9r 


+ 


uv 

r 


f  2 

9  v  +  _1_  9v 
fir2  r  8r 


(3) 


9w 

U5T  +  W 


=  -I  +  v 

9z  p  8Z 


/9Zv 

W 


8  2w  _1_  dw  +  ^  2w 


r  8r  9z 


9Zw\ 


(4) 


where  the  8/9t  and  8/9 8  terms  have  been  eliminated  by  assuming  the  flow 
is  steady  and  axisymmetric. 

The  axisymmetric  stream  function,  is  defined  so  that  Eq.  (1)  is 
everywhere  satisfied: 


u 


(5) 


In  addition,  the  dimensionless  coordinates  q  = 
and  the  quantities  O  =  ur  and  T  =  vr  are  introduced. 


(r/ rQ)2  and  % 


z/f 


( 
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Lewellen  showed  that  9r/9£  was  of  order  *  ,  where 


« 


(6) 


Therefore,  in  regions  where  the  ratio  of  stream  function  to  circulation 
is  small,  r  may  be  considered  independent  of  z,  to  order  *  .  Further, 
examination  of  Eq.  (3)  shows  that  when  v  is  independent  of  z,  the  radial 
velocity,  u,  must  also  be  independent  of  z  (unless  T  =  constant).  In  this  case, 
Eq.  (3)  reduces  to  an  ordinary,  linear  differential  equation  for  I~  in  terms  of 
u  and  v: 


r 


2  d 
37 


ur  dr 
v  37 


=  0 


(7) 


The  most  general  form  of  i)'.  consistent  with  the  condition  that  u  be 
independent  of  z,  is 

+  =  QqHUM  +  F(q)]  (8) 


so  that 


u 


w 


2Q0I 


+  F* 


) 


Lewellen^  has  further  shown  that,  since  the  equation  which  determines  ip  is 
of  order  «  ,  to  this  order  then  the  functions  f  and  F  remain  to  be  determined 
by  the  boundary  conditions. 


-7- 


It  should  be  noted  at  this  point  that  only  £  enters  in  the  determination 
of  r.  That  is,  only  that  portion  of  the  axial  flow  which  is  coupled  to  the 
radial  flow  through  the  continuity  relationship  plays  a  role. 

The  f-distribution,  or  equivalently  the  distribution  of  radial  velocity 
through  Regions  I  and  III,  will  be  determined  by  consideration  of  the  end-wall 
boundary-layer  flow.  Since  it  is  assumed  that  the  boundary-layer  mass  flow 
grows  from  zero  at  r^,  any  increase  (or  decrease)  in  the  radial  flow  in  the 
boundary  layers  must  result  in  a  corresponding  decrease  (or  increase)  in 
the  radial  flow  through  the  primary  vortex  region.  Thus,  (in  Region  I), 


f  = 


2Q 


BL 


where 


Q 


BL-  = 


ru  dz 


(H) 


(12) 


There  have  been  several  analyses  of  the  end-wall  boundary-layer 
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region.  For  the  present  study,  the  results  of  the  turbulent  analysis 

1 4 

of  Rott  will  be  used.  This  is  chosen  because,  at  the  tangential  Reynolds 
numbers  ordinarily  encountered  in  vortex  tubes,  the  boundary  layers  are 
most  likely  turbulent;  and  while  Rott' s  analysis  yields  the  simplest  expres¬ 
sion  for  the  mass  flow  in  the  boundary  layer,  the  results  agree  quite 
favorably  with  those  of  other,  more  complicated  approaches.  Assuming 
that  the  turbulent  shear  components  are  proportional  to  the  squares  of  the 
respective  velocity  components  and  that  the  coefficient  of  friction  so  defined 
is  constant,  Rott,  using  a  momentum-integral  approximation,  expressed  the 
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volume  flow  in  the  boundary  layer  in  terms  of  the  T-distribution  of  the 
external  flow  as  follows: 


where  x  =  r^  -  r,  a  =  a  velocity-profile  shape  parameter,  and  c  =  coefficient 
of  friction.  Following  Rott,  the  values 

a  =  5  ;  c  =  0.  027(Ret)“1/5  (14) 

will  be  used  for  this  study. 

The  specification  of  f  in  Region  III  is  not  quite  so  straightforward.  As 
mentioned  previously,  the  effect  of  the  discontinuity  at  the  edge  of  the  exit 
hole  (assumed  here  to  be  a  sharp-edged  orifice)  and  the  question  of 
rccirculatory  flow  along  the  axis  are  not  fully  understood.  In  addition,  there 
may  be  some  question  as  to  the  validity  of  the  boundary-layer  analysis  near 
the  axis.  All  of  these  factors  contribute  to  the  difficulty  of  specifying  f  in 
detail  in  Region  III.  Fortunately,  as  far  as  the  calculation  of  the  circulation 
is  concerned,  only  the  integral  of  f  appears  in  Eq.  (10).  Thus,  T  will  be 
insensitive  to  the  details  of  the  axial  velocity  distribution  (as  long  as  v  is 
independent  of  z).  It  will  be  assumed  then  that  a  certain  fraction,  E,,  of  the 
boundary-layer  volume  flow  at  rp  (the  radius  of  the  exit  hole),  is  ejected  and 
returned  to  the  primary  flow  to  be  redistributed  into  radial  flow  uniformly 
over  the  length  of  the  tube.  The  quantity  £  will  be  considered  parametrically 
in  the  computations,  and  its  effect  on  the  T-  and  f-distributions  will  be 
assessed.  In  addition,  inside  the  radius  of  the  exhaust,  the  simplest  possible 
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assumption  will  be  made,  namely,  that  the  radial  velocity  falls  to  zero  at 
r  =  0  linearly  with  r.  Under  these  conditions 


f  = 


(1  -  C) 


2QBL<re> 


2 

r 

re 


(15) 


when  r  <  r  . 

—  e 

Considering  Eqs  .  (11),  (13),  (14)  and  (15),  f  may  be  written 


1  >  •ne 


(16) 


f 


i  -  (i 


whe  r  e 


0  <  il  <rie 


A  =  - 


0.  27 
T7S 


(Ret) 


Zen 

<J7 


(17) 


Equations  (10)  and  (17)  determine  f  and  T  as  functions  of  the  radius  and  the 
parameters  Rer>  A,  r^/r^,  and  (»•  These  equations  have  been  solved 
iteratively  using  a  digital  computer  for  ranges  of  the  governing  parameters, 
and  the  results  are  presented  and  discussed  in  the  next  section. 
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IV.  DISCUSSION  OF  RESULTS 


The  quantity  A  characterizes  the  boundary-layer  interaction  since  it 

is  a  measure  of  the  fraction  of  the  total  mass  flow  which  passes  through  the 

boundary  layer.  If  T  were  constant  and  equal  to  Tq  for  all  r,  then  it  can  be 

seen  from  Eqs.  (13),  (14),  and  (16)  that  A  would  be  equal  to  the  fraction  of 

the  total  mass  flow  in  the  boundary  layers  at  r  =  0.  Also,  from  Eq.  (17)  it 

1/2 

is  seen  that  A  ~  !/«  ,  so  that  small  values  of  «  ordinarily  correspond  to 

large  values  of  the  boundary-layer  interaction  parameter.  Thus,  it  is 
precisely  in  the  regime  where  r  is  independent  of  z  that  the  end-wall 
boundary-layers  will  play  a  significant  role  in  the  determination  of  the  radial 
I~ -distribution. 

An  interesting  exception  occurs  in  the  limit  of  A  going  to  zero  as  a 
result  of  frictionless,  or  suitably  rotating,  end  walls  (i.  e.  ,  not  because  t 
becomes  large).  In  this  case,  the  present  problem  becomes  identical  with 

5 

that  considered  by  Einstein  and  Li.  The  stream  function  is  given  simply  by 

f  =  l  ;  n  >  ne 

f  =  ;  Ti  <  Tie  (18) 

and  Eq.  (10)  yields  the  circulation  distribution 
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where 


ne  exp(Rer/2)  -((Rer/2]+l) 

=  (Rer72)  +  1  % 


This  solution  is  shown  in  Fig.  5  for  three  values  of  the  radial  Reynolds 

number  and  r  /  r_  =  1/6. 

e  U 

Another  inte resting  limiting  case,  which  may  be  solved  analytically, 
occurs  when  Re^—  -  oo.  This  case  has  been  solved  in  Appendix  A.  In  this 
limit,  the  circulation  distribution  remains  constant  for  A  <0(1).  For 
A  >  0(1),  the  circulation  remains  constant  with  decreasing  radius  from  r^ 
until  all  of  the  flow  reaches  the  end-wall  boundary  layers.  Between  this 
radius  and  the  edge  of  the  exhaust  hole  the  total  radial  flow  remains  in  the 
boundary  layers  and  the  circulation  distribution  is  that  required  to  drive 
this  boundary-layer  flow. 

Results  of  the  digital  computer  solutions  for  T-  and  f-distributions 
corresponding  to  various  combinations  of  the  governing  parameters  are 
shown  in  Figs.  6  through  10.  All  of  these  results  were  obtained  for  a 
constant  value  of  r^/r^  =  1/6.  The  successive  figures  correspond  to 
increasing  values  of  A.  Values  of  radial  Reynolds  number  have  been  chosen 
so  that  results  representative  of  the  distributions  obtainable  are  presented. 

In  most  of  these  figures,  curves  are  plotted  only  for  4  =  0  and  4  =  1.  since 
results  for  intermediate  values  may  be  easily  inferred  from  those  presented. 
In  Fig.  7,  the  curve  for  4  =  0.  5  is  also  presented,  indicating  a  rather 
uniform  variation  of  the  distributions  with  4  between  the  extremes  of  zero  and 
one.  Results  of  the  analytical  solution  of  Appendix  A  for  Rer  =  “°°  are 
included  in  Figs.  8-10. 

The  curves  illustrate  that  when  A  <  1,  the  effect  of  the  boundary-layer 
interaction  is  relatively  small.  The  same  is  true  when  |  Re r  |  is  small 
regardless  of  the  value  of  A.  In  both  of  these  areas,  of  course,  the  influence 
of  4  is  also  small.  For  slightly  larger  values  of  A,  the  boundary-layer 
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interaction  begins  to  exert  a  considerable  influence,  especially  at  the  higher 
values  of  |Rer|.  The  effect  of  (,  also  becomes  more  pronounced  here.  A 
value  of  t,  =  1  causes  the  circulation  distribution  to  stay  closer  to  a  constant 
and  results  in  increases  of  I^/Tp  of  up  to  20  percent.  Its  influence  on  the 
mass-flow  distribution  is  more  profound  (as  seen  in  Figs.  7  and  8),  resulting 
in  negative  values  of  f,  which  implies  that  the  boundary-layer  mass  flow  is 
locally  greater  than  the  total  through-flow  so  that  radial  outflow  must  exist 
in  the  primary-flow  region.  It  is  interesting  to  note  that  in  spite  of  this 
excess  mass  flow  in  the  boundary  layers,  the  T-distribution  remains 
essentially  constant. 

At  still  larger  values  of  the  interaction  parameter,  the  F-distribution 
decreases  faster  (see  Figs.  9  and  10),  and  the  theoretical  results  for 
infinite  Reynolds  number  are  approximated  closely  by  those  for  more 
moderate  Reynolds  numbers. 

In  order  to  characterize  the  effect  of  the  interaction  in  a  simple  way, 

the  value  of  is  plotted  in  Fig.  11  as  functions  of  A  and  Re^  for 

r  /r_  =  1/6.  This  curve  also  illustrates  the  influence  of  t.  It  is  seen  that  for 
e  0 

£  =  0,  values  of  /Tq  begin  to  depart  from  one  when  A  >  1.2,  even  for  infinite 

Reynolds  number.  For  £  =  1,  this  departure  is  delayed  until  A  >  2.2.  The 

same  results  are  cros s -plotted  in  different  ways  in  Figs.  12  and  13.  Fig- 

tire  12  shows  r  /T_  vs  Re  for  constant  A,  including  the  A  =  0  results.  At 
e  0  r  & 

A  =  1  and  t,  =  1,  it  is  seen  that  T  /Tn  values  of  0.  98  are  achieved  for 

e  U 

(Re^l  >  13,  while  this  does  not  occur  for  £  =  0  until  iRe^l  >  40.  This  is  illus¬ 
trated  more  clearly  in  Fig.  13,  which  shows  curves  of  constant  as 

functions  A  and  Re  . 

v  r 

The  mass  flow  distributions  are  characterized  by  the  minimum  value 
of  f  (designated  f  )  and  the  amount  of  fluid  which  is  returned  by  the  boundary 
layers  to  the  primary  flow  (f^  -  f  ).  These  distributions  are  plotted  in 
Figs.  14  and  15.  Negative  values  of  f  indicate  reverse  flow  regions  in  the 
primary  flow. 
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In  order  to  assess  the  influence  of  rg  on  the  F-  and  f-distributions , 

solutions  were  obtained  for  specific  combinations  of  A  and  Rer  for 

re/rg  =  1/4  and  1/10,  in  addition  to  the  results  already  presented  for 

r^/r^  =  1/6.  Typical  profiles  are  shown  in  Fig.  16,  illustrating  that 

r^/Tp  and  f^  tend  to  increase  with  increasing  re/rQ  although,  at  a  given 

radius,  increasing  r  /  r-  tends  to  decrease  T.  Results  for  r  /r_  and  f 

e  0  e  0  m 

are  shown  as  functions  of  r^/r^  in  Figs.  17  and  18,  demonstrating  the 
relative  influence  of  exhaust  radius  in  various  regimes. 

In  Appendix  B,  the  effects  of  compressibility  are  considered  by 
treating  the  isothermal  case.  It  is  shown  that  the  present  results  are 
accurate  for  most  practical  applications. 


-15- 


V.  COMPARISON  WITH  EXPERIMENT  AND  COMMENTS 
REGARDING  TURBULENCE 


It  has  been  traditional  to  compare  experimental  results  and  theory  by 
matching  the  measured  circulation  distributions  (deduced  from  static 
pressure  measurements)  with  the  analytical  predictions.  In  all  cases,  it 
has  been  found  that  circulation  profiles  could  be  matched  only  when  the 
measured  laminar  radial  Reynolds  number  was  much  higher  than  its  theoret¬ 
ical  counterpart.  As  mentioned  in  the  introduction,  previous  investigations 
of  vortex  flow  in  tubes  were  based  on  one -dimensional  flow  models 
(equivalent  to  setting  A  =  0  in  the  present  theory)  and  tended  to  ascribe 
discrepancies  between  measured  and  theoretically  computed  circulation 
profiles  to  the  presence  of  turbulence  in  the  flow  field.  It  was  assumed  that 
the  turbulence  produced  a  large  "eddy  viscosity"  which  in  turn  resulted  in  a 
reduced  radial  Reynolds  number.  The  ratio  of  the  theoretically  deduced 
"effective"  radial  Reynolds  number  to  the  actual  radial  Reynolds  number 
based  on  a  laminar  viscosity  was  used  as  a  measure  of  the  turbulence  level 

Q 

in  the  tube.  Thus,  Keyes,  for  example,  was  able  to  arrive  at  a  correlation 
between  apparent  turbulence  level  and  the  tangential  peripheral  Reynolds 
number  of  the  flow. 

1  7 

Subsequently,  Kendall  and  others  pointed  out  the  significance  of  the 
end-wall  boundary-layer  interaction  and  suggested  that  the  reduction  in 
T-distribution,  which  had  been  attributed  to  turbulence,  may  indeed  be 
explainable  largely  or  wholly  by  the  three-dimensional  boundary-layer 
interaction.  It  was  argued  that  diversion  of  flow  into  the  boundary  layers 
reduced  the  radial  Reynolds  number  in  the  primary  flow,  thus  accounting  for 
the  reduced  values  of  r/Tg. 

As  the  present  study  shows,  however,  when  r/Tg  begins  to  depart 
from  one,  the  boundary-layer  growth  slows  down  and  finally  reverses, 
resulting  in  fluid  being  returned  to  the  primary  flow.  Since  only  a  fraction 
of  the  fluid's  angular  momentum  is  lost  in  the  boundary  layer  (on  the  average), 
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mass  ejection  from  the  boundary  layer  tends  to  support  the  angular 
momentum  distribution,  retarding  its  further  decline.  Finite  values  of  C 
contribute  further  to  this  trend,  resulting,  as  has  been  seen,  in  the  require¬ 
ment  that  A  be  greater  than  1.2  for  £  =  0  (or  2.2  if  £  =  1)  before  any 
significant  degradation  of  the  T-profile  occurs  due  to  the  boundary-layer 
inte  raction. 

In  order  to  determine  whether  three-dimensionality  can  indeed  explain 

all  or  part  of  the  experimentally  observed  velocity  distributions,  some  new 

experiments  were  conducted,  using  two  different  vortex  chambers.  Some  - 

1  9 

what  modified  versions  of  the  apparatus  described  by  Rosenzweig  and 

Grabowsky  and  Rosenzweig^*"1  were  used.  These  tests  will  not  be  described 

in  detail;  only  the  results  will  be  presented  here.  The  data  were  reduced 

using  the  present  theory  to  infer  a  value  of  "effective"  radial  Reynolds 

number,  after  V  and  the  value  of  A  had  been  determined  from  measured 
e 

quantities.  (Note:  Because  of  the  difficulties  associated  with  the  experi¬ 
mental  determination  of  Tq,  the  ratio  r^/Tj was  used  to  compare  with  the 
theory.)  Results  of  these  tests  are  shown  in  Fig.  19  as  a  ratio  of  actual  to 
"effective"  radial  Reynolds  numbers  vs  tangential  Reynolds  number  based  on 
tube  diameter.  A  spread  is  indicated  for  each  test  point,  the  upper  bound 
corresponding  to  £,  =  1  and  the  lower  to  £  =  0.  For  some  of  the  tests,  an 
independent  calculation  of  the  turbulent  shear  stress  is  made  possible  by 
means  of  a  torque-balance  measurement.  Such  results  can  also  be  reduced 

to  a  ratio  of  Reynolds  numbers  and  these  points  are  indicated  in  Fig.  19  as 

1  7 

well.  Also  included  are  several  hot-wire  measurements  of  Kendall.  Keyes' 

g 

leart-squares  curve  is  shown  for  comparison,  as  are  some  of  his  actual  test 
points . 

Values  of  A  in  these  experiments  ranged  from  0.6  to  1.5  and  radial 
Reynolds  numbers  varied  from  12  to  several  thousand.  Calculations 'indicate 
that  values  of  A  of  order  one  were  also  representative  of  Keyes'  experi¬ 
ments.  The  conclusion  is,  therefore,  that  early  estimates  of  turbulence 
levels  in  jet-driven  vortices,  although  based  on  one -dimensional  flow  models, 
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were  not  greatly  in  error.  Furthermore,  the  cioae  agreement  in  Reynold* 
number  ratio  between  the  result*  baaed  on  the  theory  a*  presented  herein 
and  independent  measurement*  of  turbulent  shear  lend  credence  to  the 
present  theory. 
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VI.  CONCLUSIONS 


The  equations  governing  the  interaction  between  the  primary  flow  and 
end-wall  boundary  layers  in  a  vortex  tube  have  been  presented  and  solved 
iteratively  on  a  digital  computer.  The  controlling  parameters  include  a 
boundary-layer  interaction  parameter,  A,  representing  a  measure  of  the 
mass  fraction  diverted  to  the  boundary  layers;  the  radial  Reynolds  number; 
and  a  parameter,  C, ,  representing  the  boundary-layer  mass  ejection 
occurring  as  a  result  of  the  geometrical  discontinuity  existing  at  the  exhaust- 
hole  radius.  Circulation  and  mass-flow  distributions  are  presented,  illus¬ 
trating  the  variety  of  results  obtainable  in  the  different  flow  regimes. 
Over-all  results  are  characterized  by  the  values  of  F^/Tq,  the  circulation 

ratio  at  the  exhaust  radius,  and  f  ,  the  minimum  mass  flow  in  the  primary 

m  r  7 

flow  region. 

The  results  show  that  if  A  <  1,  the  effect  of  the  boundary  layer 

interaction  on  the  circulation  distribution  is  small.  The  same  is  true  if 

|Rer  |  is  of  order  one,  regardless  of  A.  As  A  increases,  the  effect  is  more 

pronounced,  the  interaction  producing  a  reduced  r_  /  Tp  and  f ^  .  A  variation 

of  t,  from  zero  to  one  produces  a  maximum  increase  of  approximately 

20  percent  in  T^/Tp,  but  can  have  a  profound  effect  on  f^,  causing  negative 

values  (i.e.,  reverse  flow)  to  occur  in  some  instances.  It  has  been  shown 

that  in  some  cases  it  is  possible  to  have  an  excess  of  flow  in  the  boundary 

layer  (i.e.  ,  f  <  0)  while  V/V-  remains  close  to  unity. 

m  0  ’ 

The  influence  of  exhaust  radius  variation  has  been  investigated  and, 

in  addition,  analytical  solutions  for  the  limiting  cases  of  A  =  0  and 

Re  —  -  oo  have  been  presented, 
r  r 

The  theoretical  results  have  been  compared  with  experiments  by 
matching  experimental  and  theoretical  circulation  profiles  and  thereby 
determining  a  ratio  of  "actual"  (laminar)  to  "effective"  radial  Reynolds 
number.  By  comparison  with  previous  results  it  has  been  concluded  that 
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turbulence  levels  in  jet-driven  vortex  flows  are  high  and  that  estimates 
based  on  one -dimensional  assumptions  were  not  greatly  in  error. 

I 

I 


APPENDIX  A 


ANALYTIC  SOLUTION  FOR  INFINITE  REYNOLDS  NUMBER 


An  analytic  solution  to  Eqs.  (10)  and  (16)  is  possible  in  the  limit  of 
Re^-»  -  oo.  From  Eq.  (9)  it  is  seen  that  in  order  to  have  a  finite  value  of  f 
in  this  limit,  either  an  infinity  or  a  zero  in  T'  is  required.  Conversely,  any 
T-distribution  is  possible  when  f  =  0,  since  the  product  Rerf  is  then 
indeterminate.  Therefore,  it  is  evident  that  in  this  limit  T  will  remain 
constant  at  large  radii  until  all  of  the  flow  reaches  the  end-wall  boundary 
layers,  i.e.  ,  until  f  =  0.  Proceeding  inward  from  this  point,  all  the  flow 
will  remain  in  the  boundary  layers,  and  the  circulation  will  distribute  itself 
in  just  such  a  way  as  to  keep  f  =  0.  This  distribution  may  be  found  by 
writing  Eq .  (13)  in  differential  form: 


d 

37 


°BLr 


aO 


BL  dr 


dr=acrz 


(A-  1 ) 


For  Qbl  =  constant,  this  reduces  to 


1  dr  _  qc _ 

r7 37  =  qbl(i  "  a) 


( A-2) 


Equation  (A-2)  integrates  to 


r0_  A  r  , 

r  -  n  -  d)(l  -  t)  r 0  +  C2 


(A-  3) 
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in  terms  of  previously  defined  dimensionless  parameters.  The  constant  c?  is 
evaluated  by  the  condition  that  T/Tq  =  1  at  r  =  r,  the  radius  at  which  f  first 
reaches  zero.  For  r  >  r,  from  Eqs.  (11)  and  (13), 

f  =  1  -  A[1  -  (r/r0)]  (A-4) 


Therefore, 


r/rQ  =  1  -  1/A  (A-5) 

From  Eq.  (A-5)  and  the  condition  that  r/Fg  =  1,  f  =  0  at  r  =  r,  Eq.  (A-3) 
reduces  to 


r  _ _ a  -  l _ 

F0  a  -  2  +  A[l  -  (r/rQ) ] 


(A-6) 


For  4  =  0,  the  solution  for  the  T -distribution  in  the  limit  of  Re  -  oo 

A  A  T 

is  T  =  Tg  for  r  >  r  and  is  given  by  Eq.  (A-6)  for  r^  <  r  <  r.  Inside  the 

exhaust  radius  it  is  undetermined. 


The  solution  for  4  4  0  is  slightly  different.  At  the  radius  of  the 
exhaust,  flow  is  forced  out  of  the  boundary  layer  so  that  f  >  0  when  r  <  r^. 
Thus  r'  must  be  zero  here  from  Eq.  (9).  Furthermore,  differentiation  of 
Eq.  (10)  indicates  that  F*  must  remain  zero  for  r  >  r^  until 


<  0 


(A-7) 


( 
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This  in  turn  requires  f  to  be  negative  in  some  region  outside  the  radius  of 
the  exhaust.  Thus,  for  C  ^  0,  Eq.  (A-6)  applies  only  to  r  >  r  ,  with  r 
defined  so  that 


f’' 

I 

■'0 


(A-8) 


To  complete  Eq.  (A-8),  f  is  given  by 

f  =  JL[1  +  (1  .  tHfe  .1)]  ;  q  <ne 

r  __ 

f  )  :  T!  <q<q*  ( A-9) 

*o  e 


Equations  (A-8)  and  (A-9)  yield  an  algebraic  equation  relating  I~e  and  q*: 

r  r  * 

1  -  (1  -  £)[l  -  A-=^(^  -  \^)]  =  A-^t^lnl-  -  2(^  -  \GT)]  (A-  10) 

l0  l0  ■  % 

A  second  equation  relating  and  is  obtained  from  Eq.  (A-6)  which  is 
valid  at  q*: 


r 

e 


_ Q  -  1 _ 

a  -  2  +  A(1  -  Vq*) 


(A- 1 1 ) 


Equations  (A- 10)  and  (A- 11)  must  be  solved  simultaneously  to  obtain  rg  and 
q*.  A  representative  result  is  shown  in  Fig.  9. 
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APPENDIX  B 

COMPRESSIBILITY  EFFECT 


In  the  applications  and  experiments  referred  to  in  the  introduction, 
velocities  of  the  order  of  the  speed  of  sound  may  be  expected  to  occur  within 
the  vortex.  Under  such  conditions  the  assumption  of  incompressible  flow 
might  appear  to  be  a  source  of  considerable  error.  The  purpose  of  this 
appendix  is  to  prove  that  incompressible  flow  is  still  a  valid  assumption  for 
this  problem. 

A  complete  solution  of  the  compressible  problem  necessarily  involves 
a  solution  of  the  energy  equation  and  a  means  of  accounting  for  the  energy 
sources  and  sinks  within  the  flow  and  on  the  boundaries  of  the  flow.  Such 
conditions  will  vary  widely  for  the  applications  considered.  To  avoid  this 
complexity,  for  the  present  purposes  the  flow  will  be  assumed  to  be 
isothermal  (i.e.  ,  Prandtl  number  equal  to  zero).  The  isothermal  assump¬ 
tion  provides  the  most  stringent  test  of  the  effects  of  compressibility  (as 
long  as  there  are  no  heat  sources  or  sinks  present!  since  for  this  case  the 
density  change  due  to  pressure  variation  will  be  a  maximum.  For 
isothermal  flow,  it  is  also  consistent  to  keep  viscosity,  p,  a  constant.  Under 
these  conditions  the  only  one  of  Eqs  .  (1)  through  (4)  to  change  is  Eq.  (1), 
which  should  now  be  written  as 


0  0 
—  (pur)  +  —  (pwr)  =  0 
o  r  o  z 

This  calls  for  redefining  the  stream  function,  vJj,  as 

p0  ddi  P0 

pu= — ;  pw= - 

r  dz  r  or 


(B-l) 


(B-2) 
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With  this  extended  definition  of  4<>  Eqs.  (8)  and  (10)  are  unchanged  as  the 
solution  of  the  flow  in  Regions  I  and  III. 


In  Region  II,  the  assumption  of  isothermal  flow  implies  that  the 
density  gradient  across  the  boundary  layer  can  be  ignored,  while  the  radial 
variation  of  density  must  be  retained.  In  Rott' s  notation,  the  momentum 
integral  equation  in  the  tangential  direction  with  p  =  p(r)  is 


d^1 


r 


ur2(V 


-  v)  dz 


-  P^<rV>  f 


ur  dz  = 


2 

r  t 


(B-3) 


The  solution  of  Eq.  (B-3),  which  replaces  Eq.  (13),  is 


O 


BL  ” 


(B-4) 


Since  it  is  now  mass  flow  which  is  conserved  rather  than  volume  flow, 
Eq.  (11)  must  also  be  altered  to  give: 


f  = 


ZPPBL 

P0Q0* 


(B-5) 


Keeping  the  assumption  of  linear  decay  of  the  radial  mass  flow  within  the 
exhaust  radius,  and  using  Eqs.  (B-4)  and  (B-5),  Eq.  (16)  now  becomes 


f  =  1  -  A 


f  =  -3- 


(Z&i)4  [l  EW(Io)3  jn_ 

\ro  >  J  p o  W  2^ 


e 


^n>n. 


0<T,£r,e  (B-6) 


Equations  (B-6)  and  (10)  arc  the  extended  equations  replacing  Eqs. 
(16)  and  (10)  as  the  equations  to  be  solved  iteratively  on  the  computer.  The 
system  is  now  completed  by  using  Eq.  (2)  to  determine  the  pressure  and 
consequently  the  dens ity  dis t r ibution .  Neglecting  order  «,  Eq.  (2)  gives 


1  9p  _  v 

p  8r  "  r 


(B-7) 


From  the  equation  of  svate 


p  =  [p/(RT)] 


(B-8) 


Equation  (B-7)  can  be  written  as 


J_  _  _1_  r_ 

P  0r  RT  r3 


(B-9) 


which  for  isothermal  flow  integrates  to 


.2  .1 


-2-  =  exp ' 

Po 


-?/[?]' 


(B- 1 0) 


( 


with  y  the  ratio  of  specific  heats  and  Mq  the  tangential  Mach  number  at  r^. 
It  follows  directly  from  Eqs  .  (B-8)  and  (B-10)  that 


= 

Po 


(B- 11) 


The  addition  of  Eq.  (B-ll)  to  determine  p/Pq  in  Eq.  (B-6)  does  not 
greatly  complicate  the  iteration  scheme,  which  is  still  carried  through  in 
the  same  manner.  It  is  only  necessary  to  add  a  loop  determining  p/Pq  from 
the  r  obtained  in  the  previous  step.  There  is,  however,  a  new  parameter, 
YMq,  affecting  the  solutions. 

Figure  B- 1  compares  the  F  -  and  f-distributions  with  and  without  the 
compressibility  effect  for  several  sets  of  values  of  the  parameters, 
corresponding  roughly  to  the  limiting  Mach  numbers  obtained  in  the  experi¬ 
ments  discussed  in  the  paper.  Although  Mach  numbers  above  2  are  obtained, 
the  effect  of  compressibility  on  the  distributions  of  T  and  f  is  small.  It  is, 
of  course,  possible  to  increase  the  theoretical  values  of  to  the  point 
where  the  density  variation  would  cause  a  large  change  in  the  T-distribution, 
but  Mach  numbers  of  such  a  magnitude  have  never  been  obtained  experi¬ 
mentally.  Also,  an  increase  of  A,  with  constant,  will  produce  a  larger 
compressibility  effect  since  the  primary  effect  of  compressibility  is 
manifested  through  the  boundary  layer  interaction. 
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Fig.  2.  Possible  Stream¬ 
line  Patterns  in  the 
Vortex  Tube 
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3.  Vortex  Flow  Visualization  Picture  with  Dye  Injected 
into  the  Boundary  Layer,  Showing  Mass  Ejection  in- 
Vi  cinity  of  Exhaust  Radius 
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Fig.  4.  Vortex  Flow  Visualization  Picture  with  Dye  Injected 
into  the  Boundary  Layer,  with  a  Stepped-Down  End  Wall 
Showing  Mass  Ejection  in  Neighborhood  of  Step 
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Fig.  10.  Circulation  and  Stream  Function  Distribution 
as  a  Function  of  Radius  for  A  =  10 
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Fig.  13.  A  as  a  Function  of  Re  for  Constant  Circulation  at 
the  Edge  of  the  Exhaust  Hole 


Fig.  14.  Minimum  Mass  Flow  in  the  Primary  Flow 
Function  of  A  for  Constant  Re 
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Fig.  18.  Minimum  Mass  Flow  in  the  Primary  Flow 
as  a  Function  of  the  Exhaust  Radius 
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Fig.  19.  The  Ratio  of  Laminar  Radial  Reynolds  Number  to  Effective 
Turbulent  Radial  Reynolds  Number  as  a  Function 
of  Tangential  Reynolds  Number 
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Fig.  B- 1 .  Circulation  and  Stream  Function  Distributions  as  a  Function 
of  Radius  Showing  the  Effect  of  Compressibility  for  =  0.  25 
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